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Summary 

Research  conducted  under  sponsorship  of  this  grant  (F49620-98-1-0052)  focused  on 
fundamental  investigations  of  mixing,  chemical-reaction,  and  combustion  processes, 
in  turbulent,  subsonic,  and  supersonic  flows.  Specifically,  the  work  focused  on 
several  interrelated  parts: 

1.  hydrocarbon-combustion  kinetics  and  hydrocarbon-flame  behavior, 

2.  mixing  in  turbulent  transverse  jets, 

3.  mixing  in  high-speed  flows  with  a  focus  on  subsonic-diffuser  control, 

4.  Rayleigh- Taylor  instability  flows, 

5.  advances  in  Image  Correlation  Velocimetry, 
and 

6.  the  development  of  high-performance  digital  imaging  systems, 

with  significant  progress  in  all  of  these  areas,  as  documented  below. 

The  program  was  comprised  of  an  experimental  effort;  an  analytical,  modeling,  and 
numerical-simulation  effort;  and  a  diagnostics-  and  instrumentation-development  ef¬ 
fort.  Part  of  the  numerical-simulation  effort  was  in  collaboration  with  Lawrence  Liv¬ 
ermore  National  Laboratory  and  work  at  Caltech  under  ASCI/ASAP  contract  W- 
7405-ENG-48.  Parts  of  the  instrumentation-development  effort  were  cofunded  by 
AFOSR  grants  F49620-95-1-0199,  F49620-00-1-0036,  and  NSF  grants  AST9618880 
and  EIA0079871000,  in  collaboration  with  Physics,  Astronomy,  and  the  Center  for 
Advanced  Computing  Research  (CACR)  at  Caltech,  and  the  Jet  Propulsion  Labo¬ 
ratory. 
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1.  Research  description 


1.1  Hydrocarbon  flames  and  combustion  in  high-speed  flows 


1.1.1  Direct  numerical  simulations 

Detailed  numerical  investigations  on  the  ignition  of  low  molecular  weight  hydro¬ 
carbons  have  been  performed.  The  motivation  stems  from  a  proposed  concept 
( e.g .,  Edwards  1996)  for  future  hypersonic  air-breathing  combustion,  using  practi¬ 
cal  endothermic  hydrocarbon  fuels.  The  concept  pertains  to  flight  Mach  numbers 
in  the  range,  4  <  <  8.  For  typical  static  temperatures,  T0 0,  in  the  range 

250  K  <  <  300  K,  stagnation  temperatures,  To,  for  a  perfect  gas  are  expected 

to  be  in  the  range  of  1000  if  <  To  <  4000  if  (Egolfopoulos,  Dimotakis,  &  Bond 
1996,Egolfopoulos  &  Dimotakis  2000a).  Ignition  and  flame  stability  characteris¬ 
tics  of  hydrocarbons  set  the  proposed  lower  limit  of  Moo-  The  upper  limit 
is  set  by  the  cooling  capacity  of  such  fuels,  used  as  endothermic  fuels,  e.g.,  of 
kerosene-type,  to  cool  the  fuselage  and  the  internal-flow  passages.  It  has  been  es¬ 
timated  (Edwards  1996)  that  fuel  temperatures  can  be  as  high  as  1050  if,  and  for 
such  temperatures,  thermodynamics  indicate  that  fuel  molecules  are  thermally  de¬ 
composed,  through  thermal  cracking.  For  fuel  temperatures,  Tfue],  in  the  range, 
800  if  <  Tfuei  <  1050  if,  equilibrium  calculations  were  performed  (Egolfopoulos  & 
Dimotakis  1998)  to  determine  the  composition  of  the  fuel  blend  expected  from  ther¬ 
mal  cracking  of  C12H26,  a  representative  constituent  of  kerosene-type  fuels.  It  was 
found  that  the  composition  is  largely  controlled  by  nearly  equal  amounts  of  CH4 
and  C2H4,  with  smaller  amounts  of  H2,  C2H2,  C2H6,  with  C3  -  C4  hydrocarbons 
also  present. 

The  goal  of  this  investigation  was  to  numerically  assess  the  ignition  characteristics 
of  such  fuel  blends,  for  conditions  of  relevance  to  proposed  SCRAMJET  concepts. 
While  the  pressure  and  temperature  conditions  can  be  determined  with  good  accu¬ 
racy,  the  magnitude  of  strain  rates  that  will  be  locally  encountered  by  the  reacting 
layers  can  be  only  estimated  from  global  fluid  mechanics  arguments  to  be  in  the 
range  of  104  —  106s-1  (Egolfopoulos,  Dimotakis,  &  Bond  1996).  Such  strain-rate 
values  provide  guidelines  in  assessing  the  ability  of  the  fuel  blends  to  ignite. 
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Numerical  simulations  of  laminar  flame  speeds  were  performed  using  the  PREMIX 
code  (Kee  et  al.  1985).  Investigations  of  a  heated-fuel  jet  counterflowing  against  a 
heated-air  jet  permitted  an  assessment  of  the  effect  of  strain  rate.  The  non-premixed 
configuration  is  representative  of  combustion  in  SCRAMJETS,  and  its  numerical 
simulation  was  performed  using  a  stagnation-flow  code  (Egolfopoulos  &  Campbell 
1996),  modified  through  the  introduction  of  a  two-point  continuation  allowing  for 
the  description  of  the  entire  Z-curve  flame  response  to  strain  rate.  Both  codes  were 
integrated  with  the  CHEMKIN  (Kee  et  al.  1989)  and  Transport  (Kee  et  al.  1983) 
subroutine  libraries. 


Fig.  1  Effects  of  Tajr  on  Z-curve  response  of  non-premixed  C2H4/air,  for  Tfuei  = 
950 K,  using  the  WF  and  GRI  2.1  mechanisms  (see  text).  Inset:  ignition 
strain  rate  as  a  function  of  air  temperature. 


Z-curves  for  individual  constituents  and  their  blends  were  determined.  As  antici¬ 
pated,  the  ignition  propensity  of  C2H4  was  found  to  be  superior  of  that  of  CH4, 
indicating  the  need  for  further  focusing  on  ignition  promotion  of  C2H4.  The  Z-curve 
response  for  C2H4/air  is  shown  in  Fig.  1,  for  different  values  of  the  air  temperature, 
Tair.  The  large  sensitivity  of  the  ignition  strain  rate  on  Tair  is  apparent.  An  impor¬ 
tant  finding  of  this  study  was  also  that  the  omission  of  just  one  species  (vinyoxy- 
CH2CHO)  can  result  in  substantial  differences  in  the  prediction  of  ignition  strain 
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Fig.  2  Effect  of  NO  addition  on  crigri  and  <Text  of  non-premixed  C2H4/air. 

rates,  crign.  The  GRI  2.1  (Frenklach  et  al.  1995)  and  the  WF  (Wang  &  Frenklach 
1997)  mechanisms  predict  values  of  <Tjgn  that  axe  different  by  a  factor  of  30.  The 
WF  mechanism  includes  the  CH2CHO  species  while  the  GRI  2.1  does  not.  This 
study  also  demonstrated  that  the  presence  of  H2  and  C2H2  does  not  render  the 
blend  ignitable,  for  the  strain  rates  of  interest.  The  use  of  additives,  such  as  F2  and 
NO  was  found  to  increase  crign  noticeably,  as  a  result  of  the  hypergolic  nature  of 
these  two  chemicals. 

The  effect  of  additives  such  as  NO,  F2,  and  H2  on  the  ignition  of  C2H4/air,  was 
investigated,  to  assess  other  means  of  enhancing  ignition  of  C2H4.  Small  amounts 
of  NO,  F2,  and  H2  were  added  independently  into  the  fuel  and  air  streams.  It  was 
found  that  addition  of  NO  results  in  significant  increases  in  crign ,  accompanied  by  a 
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rather  moderate  decrease  of  the  extinction  strain  rate,  <jext,  when  added  to  the  air 
stream.  See  Fig.  2.  Detailed  analysis  of  the  structure  of  the  reacting  layers  showed 
that  NO  effectively  acts  as  a  catalyst  promoting  the  production  of  OH,  which  in 
turn  accelerates  the  fuel  consumption  rate  and  the  ignition  process.  The  Wang  & 
Frenklach  (1997,  WF)  kinetic  scheme  was  used  for  this  study. 

Considering  the  sensitivity  of  the  ignition  process  on  the  detailed  chemical  kinetics, 
studies  of  current  state-of-the-art  kinetics  of  C2H4  were  performed  (Egolfopoulos 
&  Dimotakis  2000b).  Six  detailed  mechanisms,  compiled  based  on  predictions  of 
experimental  results  in  homogeneous  reactors  and  flames,  were  compared.  These 
mechanisms  are:  TAN94  (Tan  et  al.  1994),  GRI21  (Frenklach  et  al.  1995),  WF97 
(Wang  &  Frenklach  1997),  MRN98  (Marinov  et  al.  1998,  Marinov  1999),GRI30 
(Bowman  et  al.  1999),  and  W99  (Wang  et  al.  1999). t 


Aerodynamic  strain  rate, a,  s1 

Fig.  3  Calculated  Z-curve  response  for  non-premixed  C2H4/air. 


The  comparative  study  included  simulations  of  both  the  Z-curve  response  and  lam¬ 
inar  flame  speeds  for  C2H4/air  mixtures.  While  laminar  flame  speed  experimental 
data  are  available,  no  ignition  and  extinction  data  are  presently  available.  Calcu- 


t  In  the  documentation  of  previous  studies,  the  acronym  WF  was  used  instead  of  WF97. 
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lated  Z-curves  are  shown  in  Fig.  3.  The  discrepancies  between  the  various  mecha¬ 
nisms  are  apparent.  In  particular,  predicted  values  of  <7ign  differ  by  factors  as  large 
as  300.  Discrepancies  of  predicted  values  for  <7ext  and  laminar  flame  speeds  are  com¬ 
parable  and  of  the  order  of  factors  of  2.3  to  3.0.  A  more-complete  documentation 
of  this  work  and  findings  can  be  found  in  Egolfopoulos  &  Dimotakis  (2000b). 

These  results  stress  the  need  for  further  kinetics  studies  of  C2H4  oxidation,  as 
well  as  the  determination  of  experimental  data  for  flame  ignition,  propagation,  and 
extinction  in  order  for  the  validations  to  be  comprehensive  and  reliable.  This  part 
of  the  work  is  a  collaboration  with  Prof.  Fokion  Egolfopoulos  of  the  University  of 
Southern  California. 

An  experimental  program  to  investigate  the  deficiencies  and  provide  data  to  address 
them  is  described  below  and  was  begun  as  part  of  this  grant  and  presently  continues 
as  part  of  AFOSR  grant  F49620-01-1-0006. 


1.1.2  Experimental  work 

The  primary  objective  of  this  part  of  the  research,  is  to  address  the  deficiencies 
in  the  experimental  data  base  and,  more  generally,  provide  a  significant  extension 
of  available  experimental  data  on  Ci-,  C2-  and  select  C3-hydrocarbon  flames,  at 
pressures  in  the  range  of  0.5  atm  <  p  <  12  atm.  Such  data  include  flame  ignition, 
propagation  and  extinction  in  laminar  premixed  flames.  Ignition  and  extinction 
strain-rates  as  well  as  laminar  flame  speeds  will  be  determined  as  a  function  of 
pressure,  temperature  and  fuel  concentration.  Comparison  between  models  and 
such  data  will  not  only  provide  a  validation  of  available  kinetic  models  but  also 
means  of  improving  them.  Selected  diagnostics  under  development  include,  but  are 
not  limited  to,  chemiluminescence  imaging,  Particle  Streak  Velocimetry  (Dimotakis, 
Debussy,  &  Koochesfahani  1981),  as  well  as  Planar  Laser  Induced  Fluorescence 
(PLIF )  for  concentration  profiles  of  radicals  such  as  CH  and  OH,  as  discussed  below. 

The  work  as  part  of  this  grant  focused  on  the  development  of  the  experimental 
diagnostics  and  configuration  at  atmospheric  pressure,  validation  of  the  technique 
against  existing  data  for  laminar-flame  properties,  and  design  and  exploratory  ex¬ 
perimentation  for  extension  to  high  pressures. 

The  experiments  entail  a  stagnation-type  jet-flow  configuration.  More  specifically, 
a  premixed  fuel-oxidizer  (air)  jet  is  generated  through  a  small-diameter  (d  ~  1  cm) 
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contoured  nozzle  that  impinges  on  a  solid  stagnation  surface.  A  planar,  nearly- 
adiabatic  flame  is  established  between  the  exit  of  the  nozzle  and  the  stagnation 
surface  at  a  location  which  is  primarily  a  function  of  the  imposed  aerodynamic 
strain-rate.  An  illustration  of  the  velocity  field  for  such  a  methane-air  flame  is 
depicted  in  Fig.  4. 


Fig.  4  Premixed  methane-air  flame  in  stagnation  flow.  Equivalence  ratio:  $  =  0.7. 
Jet-nozzle  diameter:  d  =  14mm.  Jet/plate  separation  distance:  L  =  28mm. 
Measured  jet-nozzle  exit  velocity:  f7ex  =  38  cm/s.  Hydrodynamic  strain  rate: 
cr  =  22.75  s-1. 

This  type  of  flow  has  been  extensively  documented  for  its  quasi-one-dimensional 
behavior  along  the  centerline  of  the  flow-field  (Kee  et  al.  1985).  Accurate  numerical 
descriptions  of  the  fluid  dynamics  rely  on  this  assumption.  As  a  consequence, 
significant  effort  was  expended  verifying  the  quasi-one-dimensional  nature  of  the 
flow  and  addressing  any  problems  that  make  the  flow-field  diverge  from  it.  It  was 
demonstrated  (Vagelopoulos  1999)  that  the  radial  uniformity  of  the  axial  velocity 
component  depends  on  experimental  parameters  such  as  the  design  of  the  contoured 
nozzle  and  the  ratio  of  the  separation  distance  to  the  nozzle  diameter  ( Ljd ),  as 
required  to  avoid  unfavorable  radial  pressure  gradients. 

The  jet  plenum  and  contraction  section  was  designed  to  allow  for  accurate  prediction 
of  the  boundary  layer  thickness  and  necessary  velocity  profile  corrections  at  the 
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distance  from  stagnation  plane,  cm 


Fig.  5  Sensitivity  of  flame  location  to  velocity  variation.  Methane-air  flame:  $  = 

1.0.  Jet  nozzle  diameter:  d  =  1  cm.  Jet/plate  separation  distance,  L  =  1  cm. 

Reference  jet-nozzle  exit  velocity:  Uex  =  lOOcm/s. 

the  nozzle  and  was  used  in  previous  jet-flame  work  in  our  group  (Gilbrech  1991, 
Gilbrech  &  Dimotakis  1992).  The  nozzle-exit  velocity  and  its  streamwise  gradient 
is  a  critical  boundary  condition  for  the  numerical  simulations  involved.  Figure  5 
depicts  the  sensitivity  of  the  flame  location  (represented  by  the  CH  radical  maximum 
concentration)  to  the  nozzle  exit  velocity,  as  predicted  by  numerical  simulation  (Kee 
et  al.  1985,  Bowman  et  al.  1999). 

Two  independent,  non-intrusive  methods  are  used  for  its  determination,  namely  a 
mass  flow  meter  and  a  velocity  determination  through  the  Bernoulli  pressure  drop 
in  the  jet  contraction  section.  Exploratory  experiments  verified  agreement  for  the 
two  methods  for  all  conditions  necessary  for  these  experiments,  with  an  error  of  the 
order  of,  or  less  than,  1%,  i.e.,  superior  to  typical  diagnostic  techniques  employed 
to  date.  They  were  also  consistent  with  the  measured  attendant  (small)  pressure 
rise  in  the  high-pressure  facility  (Gilbrech  1991,  Gilbrech  &  Dimotakis  1992)  that 
will  be  used  at  elevated  pressures.  Finally,  a  second  set  of  pitot-tube  measurements 
for  different  exit  velocities  verified  the  radial  dependence  of  the  nozzle-exit  velocity 
profile  to  comparable  accuracy. 

A  third  set  of  experiments  explored  the  effect  of  L/d  on  flow  uniformity,  monitored 
by  the  developed  flame  curvature  for  small  values  ( L/d  <  0.5).  As  a  general  rule, 
values  of  L/d  >  1  generated  planar  flames  (Fig.  6). 
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Fig.  6  Premixed  methane-air  flame  chemiluminescence.  Pressure:  p  —  1  atm.  Equiv¬ 
alence  ratio:  $  =  0.8.  Jet-nozzle  diameter:  d  —  1  cm.  Jet/plate  separation 
distance:  L  =  2.5  cm.  Measured  jet-nozzle  exit  velocity:  f7ex  =  75  cm/s. 
Nozzle  lip  visible  on  lower  right.  Image  resolution:  10242  pixels. 


The  stagnation  plate  was  aerodynamically  designed- to  minimize  flame  curvature 
and  instabilities.  Comparing  Fig.  4,  for  which  the  jet  flow  was  generated  through  a 
nozzle  and  a  stagnation  surface  that  were  not  optimized,  with  Fig.  6  the  substantial 
reduction  of  flame  curvature  around  the  centerline  is  evident. 

As  noted  above,  another  critical  boundary  condition  is  the  axial  velocity  gradient 
at  the  exit  of  the  nozzle.  Figure  7  depicts  the  sensitivity  of  the  flame  location  to 
this  gradient.  It  was  shown  experimentally  and  numerically  that,  for  large  L/d 
(>  1.25),  the  flow  has  jet  characteristics  in  the  vicinity  of  the  exit  of  the  nozzle, 
with  a  near-zero  velocity  gradient. 
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Fig  .  7  Sensitivity  of  flame  location  to  velocity  gradient  at  the  exit  of  the  nozzle. 
Premixed  methane-air  flame:  $  =  1.0,  p  =  latm.  Jet-nozzle  diameter: 
d  =  1  cm.  Jet-plate  separation  distance:  L  =  1  cm.  Numerical-simulation 
jet-nozzle  exit  velocity:  Uex  =  100  cm/s. 

The  axial  velocity  profile  along  the  centerline  was  determined  experimentally  through 
LDV  point  measurements  and  compared  with  two-dimensional  direct- numerical  sim¬ 
ulations  using  an  axisymmetric  code  developed  in-house  by  Ron  Henderson.  Figure 
8  depicts  the  jet-like  behavior  and  good  agreement  between  the  numerical-simulation 
and  experimental  data. 

The  lead  investigator  of  these  experiments  is  C.  Vagelopoulos,  in  collaboration 
with  J.  Bergthorson.  This  work  is  presently  in  progress  under  AFOSR  sponsor¬ 
ship  (F49620-01-1-0006)  with  studies  of  methane  and  ethylene  fuel-air  mixtures, 
before  other  hydrocarbon  blends  are  investigated. 

As  part  of  the  experimental  work  on  flames,  an  effort  was  focused  to  determine 
additional  requirements  these  experiments  will  place  on  an  imaging  system.  One 
goal  is  to  excite  CH  radicals  within  the  flame  with  a  tunable  dye  laser  operating 
in  the  UV,  at  389  nm.  CH  radicals  excited  at  389  nm  fluoresce  at  432  nm  and  this 
light  will  be  imaged. 

Such  imaging  is  difficult  because  the  fluorescence  emission  is  short-lived  and  spec¬ 
trally  amidst  the  broadband  chemiluminescence  emission  of  the  flame.  This  repre¬ 
sents  ‘noise’  that  cannot  be  removed  by  bandpass  or  other  optical  filters.  The  CH 
emission  follows  shortly  after  the  laser  excitation  and  lasts  on  the  order  of  10  ns. 
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Fig.  8  Comparison  of  experimental  (solid  symbols)  and  numerical  (open  symbols) 
data  for  cold  jet  flow.  Velocity  measured  with  LDV  diagnostics.  Numerical- 
simulation  using  axisymmetric  flow  DNS  (see  text).  Jet-nozzle  diameter: 
d  =  3  cm.  Jet/plate  separation  distance:  L  =  8  cm. 

When  using  an  ungated  CCD  imager,  chemiluminescence  is  recorded  throughout 
the  exposure.  The  short-duration  fluorescence  signal  is  negligible  when  compared 
to  this  integrated  light  from  the  flame  and  results  in  low  signal-to- noise  ratio  (SNR) 
data.  To  increase  SNR,  gating  in  time  is  required.  To  ‘shutter’  a  CCD  at  the  rel¬ 
evant  nanosecond  time  scales,  an  Intensified  CCD  (ICCD)  is  required.  Presently 
available  ICCD’s,  however,  are  characterized  by  significantly  lower  SNR’s  than  our 
in-house  digital-imaging  capability,  as  well  as  offering  lower  framing  rates  (3  fps,  or 
less).  Exploratory  experiments  were  conducted  that  confirmed  our  calculations. 

Pending  commercial  or  in-house  improvements  in  ICCD’s,  with  better  resolution 
(>  1024  x  1024  pixels)  and  framing  rate  (>  10  Hz)  diagnostics  will  primarily  rely  on 
high-resolution  chemiluminescence  imaging,  in  conjunction  with  velocity-measuring 
diagnostics.  These  provide  an  adequate  marker  for  flame  location.  At  the  high 
anticipated  pressures,  these  may  also  prove  to  be  the  preferred  option.  Gated, 
digital-imaging  design  studies  are  continuing,  under  F49620-01-1-0006,  with  an 
eye  to  improving  on  what  is  presently  available. 

This  work  is  part  of  the  Ph.D.  research  of  J.  M.  Bergthorson,  in  collaboration  with 
C.  Vagelopoulos,  and  W.  Pun,  S.  Palm,  and  Prof.  F.  Culick  at  Caltech. 
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1.2  Mixing  in  turbulent  transverse  jets 

The  mixing  behavior  of  turbulent  jets  discharging  into  a  crossflow,  or  transverse  jets, 
was  experimentally  studied  during  this  reporting  period.  Transverse  jets  are  impor¬ 
tant  to  VTOL/STOL  aerodynamics,  and  blade-and-endwall  cooling  in  gas-turbine 
engines.  The  transverse  jet  has  also  been  proposed  as  a  geometry  for  fuel  injec¬ 
tion  in  high-speed  airbreathing  propulsion  devices,  such  as  SCRAMJETS  (Mathur 
et  al.  1999).  The  measurements  of  scalar  fields  undertaken  have  enabled  analysis 
of  many  aspects  of  mixing  in  liquid-phase  (high-Schmidt-number)  transverse  jets. 
The  results  of  studies  on  Reynolds-number  effects  and  flow  dependence  of  turbulent 
mixing  are  discussed  below. 

High-resolution,  high  signal-to-noise-ratio  images  of  the  scalar  field  in  transverse 
jets  were  recorded  using  laser-induced  fluorescence  (LIF)  and  digital-imaging  tech¬ 
niques.  For  these  liquid-phase  experiments,  rhodamine-6G  was  used  as  the  flu¬ 
orescent  marker,  and  the  fluid  Schmidt  number,  Sc  =  is  estimated  to  be 

Sc  ~  8.4  x  103  (Axelrod  et  al.  1976,  Walker  1987).  With  appropriate  calibration 
and  normalization,  these  images  yield  quantitative  measurements  of  the  jet-fluid 
concentration  field  in  two-dimensional  slices.  Measurements  were  made  in  stream- 
wise  (plane-of-symmetry)  slices  of  the  jet,  as  shown  in  Fig.  9.  Images  were  also 
recorded  in  a  plane  normal  to  the  freestream  flow,  as  seen  in  Fig.  10.  These  images 
were  recorded  in  the  far-field  of  this  flow,  at  x/dj  =  50.  The  Reynolds  numbers 
investigated  were  1.0  x  103  <  Rej  <  20  x  103,  and  the  jet-to-freestream  velocity 
ratio,  VT  =  Uj/Uoo,  was  Vr  =  10. 

The  statistics  of  turbulent  mixing  in  the  transverse  jet  are  described  in  terms  of 
the  probability-density  function  (PDF)  of  the  concentration  field.  PDFs  define  the 
statistics  of  fluctuating  quantities  (scalar  concentration,  in  this  case),  and  their 
moments,  and  axe  the  basis  for  statistical  approaches  toward  turbulent  flow  with 
chemical  reactions  (Pope  1985,  Goldin  &;  Menon  1997).  For  the  transverse  jet, 
spatial  PDFs  of  jet-fluid  concentration  are  computed  from  a  sequence  of  512  (1024  x 
1024)-pixel  images,  yielding  ~  5  x  108  individual  (point)  measurements  of  the  scalar 
(jet-fluid  concentration)  field.  Jet-fluid  concentration  PDFs  in  the  far-field  (x/dj  = 
50)  of  the  jet  are  shown  in  Fig.  11  for  Rej  =  1.0, 2.0, 5.0, 10,  and  20  x  103.  The  PDFs 
are  monotonic-decreasing  functions  at  low  Reynolds  numbers,  but  show  a  growing 
peak  with  increasing  Reynolds  number.  At  the  higher  Reynolds  numbers,  there 
is  decreased  probability  for  unmixed,  high-concentration  fluid,  and  a  developing 
preference  for  a  particular  mixed-fluid  concentration,  Cpeak-  These  experiments 
and  analyses  indicate  more-uniform  mixing  ( i.e .,  better  spatial  homogenization) 
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FlG .  9  Laser-induced  fluorescence  image  data  in  a  liquid-phase  transverse  jet .  Stream- 
wise  (plane-of-symmetry)  slice  of  jet  for  VT  =  10.  Image  intensity  scaled  by 
x1/2  (for  display  purposes  only)  to  compensate  for  downstream  decay.  Left: 
Re j  =  1.0  x  103.  Right:  Re j  =  10  x  103. 


FlG.  10  Jet-fluid  concentration  in  a  transverse  slice  of  jet  at  downstream  location 
x/dj  =  50  for  VT  ~  10.1.  Left:  Rej  =  1.0  x  103.  Right:  Rej  =  10  x  103 


with  increasing  Reynolds  number,  in  the  range  1.0  x  103  <  Re j  <  20  x  103. 

Comparison  of  the  transverse-jet  PDF  with  the  PDF  of  ordinary  turbulent  jets  re¬ 
veals  flow-dependent  differences  in  turbulent  mixing.  Previous  work  in  the  group 
focused  on  the  mixing  behavior  of  liquid-phase,  axisymmetric,  turbulent  jets  dis¬ 
charging  into  a  quiescent  reservoir  (Catrakis  &  Dimotakis  1996).  The  PDFs  of  these 
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Fig.  11  Transverse-jet  PDFs  at  x/dj  =  50.  Concentrations  normalized  by  plenum- 
exit  concentration,  C  —  c/cq.  Increasingly  peaked  PDFs  at  C  ~  0.04  with 
increasing  Reynolds  number  for  Rej  =  1.0, 2.0, 5.0, 10,  and  20  x  103. 

turbulent  jets,  plotted  in  Fig.  12,  show  decreasing  peaks  with  increasing  Reynolds 
number:  the  “valleys”  fill  in,  and  the  PDFs  tend  toward  monotonic-decreasing 
functions  with  increasing  Reynolds  number.  This  is  substantially  different  from 
the  PDFs  of  transverse  jets,  which  show  growing  peaks  with  increasing  Reynolds 
number.  Turbulent  mixing  is  found  to  be  non-universal  (flow-dependent),  based  on 
the  difference  between  transverse  jets,  and  ordinary  turbulent  jets. 

Two-dimensional  image  measurements  of  the  scalar  field,  recorded  in  LIF  images 
like  those  shown  in  Fig.  10,  enable  field  measurements  of  scalar  differences.  Scalar 
differences, 

Ar(7  =  C(x  +  r,t)  -  C(x,t)  ,  (1) 

also  known  as  scalar  increments,  are  differences  between  simultaneous  measurements 
of  the  scalar  field  at  two  points  separated  by  a  distance  r.  The  statistics  of  scalar 
increments  are  important  in  turbulence  theory  ( e.g .,  scalar  structure  functions), 
and  have  been  studied  with  hot-wire  measurements  in  a  variety  of  flows  (Antonia 
et  al.  1984,  Ching  1991,  Mydlarski  &  Warhaft  1998).  Field  measurements  (rather 
than  point  measurements)  of  scalar  increments  are  made  by  spatially  shifting  the 
measured  concentration  fields,  and  subtracting  from  unshifted  images.  Examples 
of  two-dimensional  scalar  difference  fields  for  horizontal  shifts  of  2  pixels,  and  16 
pixels,  are  shown  in  Fig.  13. 


Fig.  12  Jet- in-quiescent-reservoir  PDFs.  Lines  of  increasing  solidity  denoting  in¬ 
creasing  Re j  =  4.5  x  103,  9  x  103,  and  18  x  103  (Catrakis  &  Dimotakis 
1996,  Fig.  8). 


Fig.  13  Two-dimensional  scalar-difference  field.  Left:  Horizontal  increment  of  2 
pixels.  Right:  Horizontal  increment  of  16  pixels. 


PDFs  of  scalar  increments  are  computed  from  normalized  histograms  of  the  scalar- 
difference  images.  The  PDFs  of  scalar  increments  at  x/dj  =  50  for  a  variety  of 
separation  distances  are  shown  in  Fig.  14.  At  a  jet  Reynolds  number  of  Re}  = 
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Fig.  14  PDF  of  scalar  increments  at  x/d-}  =  50  for  horizontal  separations  of  2,  4, 
16,  64,  and  128  pixels.  Innermost,  triangular  PDF  is  for  smallest  separa¬ 
tion  distance,  while  outermost,  broad  PDF  is  for  largest  separation.  For 
comparison,  the  scalar  diffusion  scales  for  the  Re }  —  2.0  x  103  and  20  x  103 
jets  are  3  and  0.5  pixels,  respectively.  Top:  Rej  =  2.0  x  103.  Bottom: 
Re j  =  20  x  103. 

20  x  103,  the  PDF  shows  exponential  tails  for  small  separations  (r  =  2  and  4  pixels). 
As  the  separation  distances  increase,  the  tails  of  the  PDFs  decay  faster  than  an 
exponential.  Others  have  reported  Gaussian  PDFs  for  large  separation  distances 
( e.g .,  Ching  1991).  At  lower  Reynolds  numbers,  Rej  =  2.0  x  103,  the  PDFs  are 
wider,  indicating  less  mixing  (poorer  spatial  homogenization).  The  PDFs  for  small 
separations  are  stretched  exponential,  and  the  tails  of  the  PDFs  once  again  decay 
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faster  than  an  exponential  for  large  separation  distance.  These  results  are  consistent 
with  those  of  investigations  in  grid  turbulence,  Rayleigh-Benard  convection,  and 
other  flows. 

This  work  is  part  of  the  Ph.D.  research  of  Jerry  Shan.  A  presentation  of  these 
results  was  made  at  the  recent  APS  meeting  of  the  Division  of  Fluid  Dynamics 
(Shan  &  Dimotakis  2000).  A  more  complete  documentation  is  available  in  Shan 
(2001)  and  Shan  &  Dimotakis  (2001). 


1.3  Mixing  in  high-speed  flows 

High-speed  internal  flows  are  particularly  sensitive  to  flow  geometry  as  well  as  skin- 
friction  effects.  Traditional  approaches  to  addressing  the  difficulties  that  ensue 
rely  on  variable  geometry,  with  high  attendant  penalties  in  weight  and  complexity. 
The  part  of  the  work  addressing  high-speed  mixing  focused  on  applications  of  a 
flow  field,  similar  to  that  anticipated  in  supersonic/hypersonic  propulsors,  shown  in 
Figs.  15  and  16.  During  the  course  of  this  grant,  the  subsonic  behavior  of  such  flows 
was  investigated.  This  study  has  important  implications  in  the  control  of  diffuser 
performance  through  aerodynamic  means. 


Fig.  15  Schematic  of  supersonic  expansion  turn  flow,  with  no  mass  injection  from 
the  porous  wall. 
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Fig.  16  Schematic  of  supersonic  expansion  turn  flow,  with  mass  injection  from  the 
porous  wall.  Boundary  layer  is  now  lifted  from  the  wall  to  form  a  shear 
layer,  entraining  the  transpiration  (lower)  injection  stream. 

As  opposed  to  supersonic  flow,  where  the  flow  will  remain  attached  over  an  ex¬ 
pansion  ramp,  in  subsonic  flow,  the  flow  separates  forming  a  recirculation  zone 
downstream  of  the  ramp.  The  control  of  this  recirculation  zone  through  mass  injec¬ 
tion  was  investigated,  taking  advantage  of  variable-geometry  effects  with  no  moving 
parts. 


SuDcc''! 


Fig.  17  Supersonic  expansion-turn  plate  assembly. 
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Fig.  18  Supersonic  expansion-turn  perforated  plate  (ramp). 

A  new  splitter  plate  and  lower  section  nozzle  block  were  designed  and  fabricated  as 
part  of  this  project  for  use  in  the  GALCIT  Supersonic  Shear  Layer  (S3L)  facility. 
Figure  17  shows  the  total  assembly  of  the  splitter  plate,  nozzle  block,  perforated 
plate,  and  lower  guidewall.  A  perforated  plate  was  designed  for  use  as  the  permeable 
inclined  wall  in  the  schematics  depicted  in  Figs.  16  and  17  and  is  shown  in  Fig.  18. 
It  has  a  regular  array  of  3611  holes  of  diameter  0.065  in,  with  an  open-area  ratio 
approximately  65%,  chosen  to  avoid  the  jet-coalescence  instability  documented  by 
Loehrke  &  Nagib  (1972). 

A  first  set  of  experiments  was  conducted  with  the  perforated  plate  replaced  by  a  solid 
ramp  to  determine  the  reference  behavior  of  the  flow  with  no  mass  injection.  Figure 
19  plots  normalized  streamwise  velocity  profiles  for  different  maximum  freestream 
velocities.  Figure  20  shows  a  sketch  of  the  variation  in  the  velocity  profile  as  the 
upstream  velocity  is  increased.  Figure  21  plots  the  normalized  velocity,  at  the  fixed 
measuring  station,  versus  upstream  Reynolds  number,  for  the  three  lowest  pitot 
pressure  probes  in  the  test  section.  This  plot  suggests  that  the  location  of  the 
reattachment  point  for  this  flow  varies  only  weakly  with  Reynolds  number.  This  is 
not  in  accord  with  the  conclusion  in  Song,  DeGraaff,  &  Eaton  (2000)  who  report  a 
reattachment  point  that  is  independent  of  Reynolds  number. 


A  second  set  of  experiments  was  then  performed  with  injection  of  a  density-matched 
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FlC.  19  Normalized,  streamwise  velocity  profile  for  one-stream  flow  (solid  ramp). 


Station  2 
(pitot  tubes  and 


Fig.  20  Schematic  of  subsonic  flow  field  with  no  mass  injection,  showing  recircula¬ 
tion  zone  and  downstream  measuring  station 


helium/argon  mixture  into  an  upper,  high-speed  stream  of  nitrogen.  Figure  22 
plots  the  normalized  velocity  profile  for  various  levels  of  mass  injection  of  a  density- 
matched  He/Ar  mixture  into  the  upper,  N2  stream.  A  run  was  performed  with 
no  mass  injection  (solid  inclined  surface),  and  compared  to  a  run  with  the  same 
upstream  velocity,  but  with  the  small  mass  injection  derived  from  the  initial  volume 
of  gas  in  the  low-speed  plenum  (initially  at  a  slightly  higher  pressure) .  Subsequent 


Fig.  21  Normalized  velocity  on  three  lowest  pitot-pressure  measuring  probes  at  the 
fixed  measuring  station,  vs.  Re. 


Fig.  22  Normalized  streamwise  velocity  profiles  for  mass  injection  of  a  density- 
matched  He/Ar  mixture  (see  text). 

runs  were  performed  with  mass  injection  from  the  low-speed  stream  at  velocities  of 
1.8,  6.0,  and  18  m/s,  respectively.  In  all  cases,  the  velocity  of  the  upper  stream  was 
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Fig.  23  Schlieren  visualization  of  mixing  after  expansion  turn;  U\  ~  55  m/s  [N2] 
and  U2  —  1.8 m/s  [He/Ar]. 

held  constant,  the  value  in  the  legend  representing  the  maximum  velocity  achieved 
at  the  downstream  measuring  station.  As  expected  and  inferred  from  these  data, 
the  reattachment  point  moves  downstream  as  the  mass  injection  rate  is  increased. 


Fig.  24  Schlieren  visualization  of  mixing  after  expansion  turn:  U\  ~  55  m/s  fNol, 
U2  ca  6 m/s  [He/Ar], 


Schlieren  visualizations  were  recorded  with  a  10242-pixel,  30fps  CCD  camera  by 
Silicon  Mountain  Design,  Model  SMD  1M30.  Figures  23  and  24  show  the  upstream 
flow  with  U\  ~  55  m/s,  and  U2  —  1.8  m/s  and  6  m/s,  respectively.  In  these  fig¬ 
ures,  a  small  region  of  uniformly  mixed  fluid  is  discernible,  just  downstream  of 
the  perforated  plate,  with  a  recirculation  zone  downstream  of  this  region.  Figure 
25  is  a  composite  of  two  runs  performed  with  the  same  inflow,  U\  ~  55  m/s  and 
U2  ~  18  m/s,  one  with  the  upstream  region,  the  other  with  the  downstream  region. 
In  this  figure,  the  transition  from  a  separated  region  towards  a  classical  shear  layer 
is  apparent,  with  a  separated-flow  reattachment  beyond  the  measuring  station. 


Fig.  25  Schlieren  visualization  of  mixing  after  expansion  turn:  U\  ~  55  m/s  [N2], 
U2  —  18  m/s  [He/Ar]. 


Many  aspects  of  diffuser  operation  can  be  parameterized  in  terms  of  a  pressure 
coefficient, 

r  _  P2  -Pi 

p  ipu*  ' 

This  is  plotted  in  Fig.  26,  with  Station  1  located  on  the  top  guidewall,  above  the  start 
of  the  ramp,  and  Station  2,  also  on  the  top  guidewall,  at  the  downstream  location 
of  the  total-pressure-probe  array,  as  indicated  in  Fig.  20.  The  two  data  points  at 
U2  =  0  correspond  to  flow  with  no  mass  injection  and  small  mass  injection,  from 
gas  already  in  the  lower  plenum,  whose  pressure  falls  in  the  course  of  a  run,  as 
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described  above.  As  expected,  the  pressure  coefficient  decreases  with  increasing 
mass  injection,  as  the  reattachment  point  is  pushed  further  downstream  and  the 
flow  evolves  to  form  a  classical  shear  layer  with  a  small  streamwise  pressure  gradient 
( e.g .,  Dimotakis  1991). 

This  work  is  part  of  the  graduate  research  by  Wei- Jen  Su  and  Michael  Johnson. 
It  is  a  precursor  to  research  targeting  the  performance  and  control  of  supersonic 
internal  flows,  as  would  be  encountered  in  hypersonic  propulsors.  This  is  presently 
in  progress,  as  part  of  work  undertaken  under  AFOSR  Grant  F49620-01-1-0006. 


1.4  Rayleigh- Taylor  instability  flow 

The  Rayleigh- Taylor  instability  (RTI)  occurs  whenever  variable-density  flows  are 
subjected  to  acceleration  in  a  direction  opposite  that  of  the  density  gradient,  e.g., 
as  occurs  when  a  heavy  fluid  finds  itself  atop  a  lower-density  fluid  in  a  downward 
gravitational  acceleration  (e.g.,  Rayleigh  1883;  Taylor  1950;  Chandrasekhar  1955, 
1961;  Duff  et  al.  1962,  Sharp  1984).  RTI  flows  arise  in  a  variety  of  contexts,  such  as 
in  combustion,  rotating-machinery  flows,  in  inertial-confinement  fusion,  supernovae- 
explosion  flows,  geophysical  flows,  and  others.  The  flow  eventually  becomes  turbu¬ 
lent  as  secondary  (Kelvin-Helmholtz)  instabilities  develop,  increasing,  in  time,  the 
spectrum  of  spatial  and  temporal  scales  that  participate  in  the  dynamics. 

Simulation  and  modeling  of  RTI  flows  represents  an  important  test  case  for  hy¬ 
drodynamic  codes,  exercising  many  aspects  of  turbulence,  transport,  and  diffusion 
in  nonuniform-density  flows  subject  to  external  body  forces.  As  opposed  to  ho¬ 
mogeneous  turbulence  in  uniform-density  flows,  RTI  flows  are  driven  by  a  directed 
internal-force  term  that  is  capable  of  sustaining  anisotropy  in  fully  developed  flow. 
Diffusive  mixing,  which  anneals  composition  (density)  differences,  plays  a  dynamic 
role  in  this  flow  and  must  be  captured  to  represent  the  dynamics. 

RTI  models  typically  assume  that  the  extent,  h,  of  the  turbulent  mixing  zone  in  the 
direction  of  the  acceleration,  g,  following  an  initial,  linear-instability,  exponential- 
growth  period,  grows  quadratically  in  time  (for  constant  g  =  |g|),  i.e., 

h(t)  ~  aAgt2  ,  (3) 

where  a  is  accepted  as  an  empirical  constant  (0.02  <  a  <  0.08),  and 

A  —  P2  ~  Pi  _  P2/P1  ~  1 
f>2  +  Pi  P2/ Pi  +  1 


(4) 


26 


is  the  Atwood  number,  with  p\  and  pi  the  densities  of  the  light  and  heavy  fluids 
across  the  interface,  respectively. 


Other  parameters  for  this  flow  are  the  Schmidt  number  (ratio  of  the  kinematic 
viscosity,  v ,  to  the  species  diflusivity,  V)  and  the  Reynolds  number,  i.e., 


with  V  the  species  diflusivity,  and 


(5) 


„  pLU  LU 

Re  =  -  =  -  , 

fx  v 

the  (outer)  Reynolds  number,  which,  for  this  flow,  may  be  defined  as, 


(6a) 


„  hh 
Re  —  — 


(6b) 


A  host  of  experimental  evidence  as  well  as  theoretical  arguments  indicate  that  it 
must  be  at  least  1  —  2  x  104,  or  so,  for  the  flow  to  be  regarded  as  bona  fide  turbulent 
(Dimotakis  1993,  2000). 


Direct  numerical  simulation  of  RTI  flows  were  performed  under  the  sponsorship  of 
this  AFOSR  grant,  as  a  collaborative  effort  with  the  Caltech  ASCI/ASAP  program 
and  the  Lawrence  Livermore  National  Laboratory  (LLNL).  A  pilot  experimental 
effort  was  also  undertaken,  as  will  be  described  below. 


The  numerical  simulation  effort  focused  on  a  study  of  the  incompressible  3-D  RTI 
flow  field,  computed  on  the  multiprocessor  ASCI  machines  at  LLNL  by  A.  Cook 
(LLNL),  followed  by  quantitative  data  analysis  and  3-D  computer  visualization  by 
S.  Deusch,  S.  Lombeyda,  J.  Patton,  and  P.  Dimotakis,  at  Caltech.  Two  sets  of  sim¬ 
ulation  runs  were  performed,  with  the  second  set  attaining  a  maximum  Reynolds 
number  of  «  3700,  i.e.,  roughly  a  third  of  the  mixing-transition  value  of  Re  «  104 
(Dimotakis  1993,  2000).  An  illustration  of  the  resulting  geometry  of  the  intermedi¬ 
ate  isosurface  from  the  latter  simulations  is  depicted  in  Fig.  27. 

The  boundary  conditions  for  these  simulations  were  space-periodic  on  the  vertical 
faces  of  the  computational  domain,  with  no-slip  on  the  top/bottom  walls.  Three 
cases  were  run,  with  different  initial-perturbation  spectra  with  peak  wavenumbers 
lower  than  the  most-unstable  mode  according  to  the  linear-stability  analysis  of  Duff 
et  al.  (1962),  for  Case  A,  comparable,  for  Case  B,  and  higher  for  Case  C.  The  initial 
perturbation  spectra  were  normalized  to  the  same  mean-square  value  for  all  three 


runs. 
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Fig.  27  Time-evolution  of  intermediate,  (p  =  2)-isosurface  in  Rayleigh- Taylor  flow. 
Pure  heavy  (p  =  3)  fluid  is  red,  pure  light  (p  =  1)  fluid  is  blue,  and  the 
(p  =  2)-isosurface  is  green.  Times  for  the  three  images  displayed:  t/r  —  0.0, 
3.44,  and  4.63,  with  r  =  y/L/Ag  and  L  the  horizontal  extent  of  the  flow 
domain  (Cook  &  Dimotakis  2000a). 


Figure  28  depicts  the  overall  extent,  h(t),  of  the  RTI  mixing  zone,  along  with 
indicatrices  corresponding  to  t}’  ~  and  t~  growth.  Growth  during  the  diffusion- 
dominated,  h  oc  it1/2,  regime  is  the  same  for  all  cases,  i.e...  independent  of  the 
initial  perturbations,  and  in  good  agreement  with  analytical  predictions  for  purely 
diffusive  growth.  The  onset  of  a  second  growth  stage  is  conspicuous,  with  growth 
that  is  faster  (higher  power  of  t)  than  diffusive.  This  occurs  at  a  different  time  for 
each  of  the  three  cases,  reflecting  differences  in  the  initial  modal  seeding  in  each 
case.  Available  linear-stability  analysis  (Duff  el  al.  1962)  does  not  account  for  the 
order  in  which  the  growth  breaks  out  of  the  diffusive  regime.  Also,  as  can  be  seen, 
only  Case  C,  which  was  seeded  with  the  highest  wavenumber  initial  disturbances, 
appears  to  grow  approximately  quadratically.  Cases  A  and  B  do  not  support  the 
conventional  modeling  of  this  phenomenon. 

Analysis  of  the  results  yielded  statistics  for  Taylor  microscales,  indicating  that  the 
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Fig.  28  Mixing-zone  extent,  h/L,  vs.  t/r,  for  the  three  cases  (Cook  &  Dimotakis 
2000a). 


directed  buoyancy  force  that  drives  this  flow  sustains  an  anisotropic  field,  precluding 
the  use  of  standard  turbulence  theory.  Values  were  found  to  be  similar  in  scaling 
and  magnitude  to  ones  measured  in  high  Reynolds  number  turbulent  jets. 

A  summary  presentation  of  results  from  this  work  was  made  at  the  ICTAM  2000 
meeting  (Cook  &  Dimotakis  2000b).  A  first  documentation  of  this  work  was  issued 
as  an  LLNL  report  and,  revised,  subsequently  submitted  for  publication  to  the  J. 
Fluid  Mechanics  (Cook  &  Dimotakis  2000a). 

A  preliminary  experimental  investigation  of  the  Rayleigh- Taylor  instability  was  also 
undertaken  as  part  of  this  work.  The  investigation  spans  the  experimental  process 
from  conceptual  and  experimental  design  to  the  execution  of  a  set  of  experiments 
tailored  to  capture  a  flow  that  is  very  difficult  to  probe. 

The  design  and  fabrication  of  a  new  pilot  facility  was  completed.  This  design 
adapted  and  improved  upon  an  established  method  for  keeping  two  fluids  in  an 
unstable  configuration  prior  to  flow  initiation  (Dalziel  et  al.  1999).  Two  experiments 
were  successful,  at  an  Atwood  number,  A  —  1.26  x  10-4  (Eq.  4),  with  results  in 
reasonable  accord  with  previous  experiments. 
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Fig.  29  Laser-induced  fluorescence  images  of  Rayleigh-Taylor  instability  (. A  -  1.26  x 
10~4).  Experiment  A. 

Flow  visualization  was  performed  using  laser-induced  fluorescence  with  sodium 
fluorescein,  at  a  dye  concentration  of  4.18  x  10-4g/L  (molar  concentration  of 
1.11  x  10-6  M),  in  the  top  chamber,  and  an  argon-ion  laser  as  the  excitation  source. 
Images  were  acquired  using  the  Caltech-developed  Cassini  digital-imaging  system 
(see  Sec.  1.6). 


Fig.  30  Laser-induced  fluorescence  images  of  Rayleigh-Taylor  instability  (. A  -  1.26  x 

10~4).  Experiment  B  (see  text). 


Figure  29  shows  the  evolution  of  a  well-developed  Rayleigh-Taylor  mixing  zone. 
This  was  recorded  approximately  18  s  after  the  retraction  of  the  partition.  As  can 
be  seen,  the  instability  thickness  above  and  below  the  interface  is  approximately  the 
same,  with  the  four,  or  five  structures  arranged  in  an  alternating  pattern.  Similar 
inferences  can  also  be  drawn  from  Fig.  30. ^  Following  Sharp  (1984),  the  flow  can 


t  The  images  in  Fig.  30  were  extracted  from  a  sequence  recorded  with  the  optical  axis  below  the 
original  interface  level,  to  avoid  obtrusion  by  the  groves  that  guide  the  separation-plate  assem¬ 
bly^  c/.  Fig.  29),  while  maintaining  an  optical  axis  perpendicular  to  the  test-section  vertical 
wall,  to  avoid  perspective  distortions.  As  a  consequence,  the  developing  Rayleigh-Taylor  layer 
is  seen  near  the  top  of  the  field  of  view. 
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be  characterized  as  having  evolved  well  into  Stage  2  (penetration)  of  the  instability. 
Also  apparent  in  most  of  the  structures  are  roll-ups,  or  “mushroom  caps”,  with 
enhanced  mixing  at  the  roll-up  interfaces.  Time-histories  of  mixing-zone  height, 
top- /bottom-layer  penetration,  and  mean  thickness  of  the  developing  flow  are  shown 
in  Fig.  31.  In  12.7  s  from  the  instant  the  separator  trailing  edge  passed  the  laser 
sheet,  the  interface  height  dropped  10  mm.  The  growth  of  the  RTI  mixing  zones 
was  found  to  exhibit  an  approximately  quadratic  dependence  on  time,  in  reasonable 
accord  with  Eq.  3. 


Time,  sec 


©  Interface 

•  Top  ofthe  RT  Front 
A  Bottom  of  the  RT  Front 
o  Top  /  Bottom  Average 


Fig.  31  Displacement  time-histories  for  the  interface,  both  sides  of  the  RTI  front, 
and  their  average  ( A  =  1.26  x  10-4). 


Using  the  model  in  Eq.  3,  curve  fits  with  an  adjustable  origin  in  time,  as  required 
by  the  data,  are  shown  in  Fig.  31.  The  data  are  consistent  with  an  approximate 
value  of  the  model  constant  of  a  ~  0.6,  with  varying  values  for  the  time  origin,  to- 
This  value  of  a  is  in  reasonable  accord  with  previously  published  values. 

The  maximum  Reynolds  number  attained  in  these  experiments  (Eq.  6b)  was  modest, 
in  particular,  considerably  smaller  than  values  required  to  place  the  flow  in  the  fully 
developed  turbulence  regime  (Dimotakis  2000).  Higher  values  would  have  to  be 
attained  if  the  influence  of  turbulent  mixing  on  this  density-difference-driven  flow 
are  to  be  documented. 

This  part  of  the  work  was  performed  by  Jimmy  Fung  and  Wei-Jen  Su  and  was  partly 
supported  by  the  Caltech  ASCI/ASAP  program  (DOE  Grant  W-7405-ENG-48). 
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1.5  Advances  in  Image  Correlation  Velocimetry 

The  work  also  advanced  Image  Correlation  Velocimetry  (ICV:  Tokumaru  k  Dimo- 
takis  1995,  Gornowicz  1997),  which  extracts  velocity-field  information  from  succes¬ 
sive  images  of  any  Lagrangian  markers  (particles,  dye,  etc.).  In  a  collaborative  effort 
with  R.  Pitz  (Vanderbilt),  it  was  applied  to  Lagrangian  flow-tagging  techniques 
(Ozone-tagging  Velocimetry  (OTV),  or  hydroxyl-tagging  velocimetry  (HTV))  for 
low-  and  high-temperature  gas-phase  flows  (Pitz  et  al.  1998,  Wehrmeyer  et  al.  1999). 

In  these  techniques,  a  laser-beam  grid  is  applied  using  a  193  nm  ArF  excimer  laser 
to  produce  lines  of  ozone  in  an  air  stream  (OTV),  or  hydroxyl  radicals  in  a  burning 
hydrocarbon  flame  (HTV).  These  species  fluoresce  and  can  be  used  to  mark  advec- 
tion  of  the  grid  by  laser  induced  fluorescence  (LIF),  using  a  248  nm  KrF  laser.  ICV 
is  then  used  to  infer  the  velocity  field  from  two  realizations  of  the  advected  LIF  grid 
(Fig.  32). 


Fig.  32  Image  Correlation  Velocimetry  application  to  03-tagged  grid.  Arrows  de¬ 
note  in-plane  velocity  field  vectors  of  gas-phase  flow,  w  ~  5  —  7  m/s.  Under¬ 
lying  grey-scale  image  depicts  the  superposition  of  the  initial  and  advected 
grid  of  the  fluorescing  O3  streaks. 


In  other  related  work  in  ICV,  a  collaboration  with  6.  Savas  (Berkeley)  was  under¬ 
taken  to  provide  a  quantitative  assessment  of  ICV,  Lagrangian  Particle  Tracking 
(LPT;  Sholl  k  Savas  1997),  and  other  PIV  methods.  A  paper  documenting  this 
study  was  recently  submitted  for  publication  (Deusch  et  al.  2000).  The  work  on 
ICV  is  part  of  the  research  effort  by  S.  B.  Deusch,  in  collaboration  with  P.  Dimo- 
takis. 
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1.6  Advances  in  digital-imaging  and  high-speed  data  acquisition 

Phase  1  of  a  multi-year  project  on  high-performance  digital-imaging  development 
was  completed.  The  new  system  has  been  dubbed  the  KFS  system,  after  its  kilo- 
frame-per-second  framing-rate  specification,  and  represents  a  substantial  enhance¬ 
ment  of  our  previous-generation  Cassini  system. 

The  Cassini  system  relies  on  CCD’s  developed  for  the  NASA  Cassini  mission,  devel¬ 
oped  by  J.  Janesick,11  M.  Wadsworth,  S.  A.  Collins,  S.  T.  Elliott,  and  others,  of  the 
JPL  Space  Instrumentation  Implementation  Section.  The  Cassini  system  camera 
head  and  data  acquisition  electronics  were  designed  and  integrated  by  Dan  Lang, 
in  collaboration  with  Pavel  Svitek  and  Paul  Dimotakis  at  Caltech.  The  Cassini 
system  is  capable  of  lOfps  (frames  per  second),  at  up  to  12  bits/pixel  at  10242- 
pixel  frame  resolution,  or  20  fps  at  5122  -pixel  frames,  which  even  though  completed 
several  years  ago,  remains  ahead  of  commercially  available  systems,  in  terms  of  its 
combination  of  spatial  resolution,  dynamic  range  (SNR),  and  framing  rate.  It  was 
used  in  the  transverse-jet  mixing  studies  (Sec.  1.2)  and  is  presently  in  use  in  the 
hydrocarbon-flame  studies  (Sec.  1.1.2),  under  sponsorship  of  the  current  AFOSR 
F49620-01-1-0006  grant. 

Initially  under  AFOSR/DURIP  Grant  F49620-95-1-0199,  Phase  1  of  the  KFS  sys¬ 
tem  was  completed  as  part  of  this  program,  with  some  additional  support  from  the 
AFOSR  Grant  F49620-00-1-0036  on  Aerooptics  studies,  with  which  the  develop¬ 
ment  and  use  of  this  equipment  is  shared.  It  should  be  noted  that  digital  imaging 
has  been  part  of  all  the  experimental  efforts  under  this  grant.  Those  continuing 
under  current  AFOSR  sponsorship  (Sec.  1.2)  await  the  next  developments  for  their 
completion. 

The  KFS  system  is  designed  to  record  10242-pixel  frames,  at  12-bits/pixel,  at  rates 
as  high  as  103  fps  (frames/sec),  for  an  aggregate  data  rate  as  high  as  109  measure¬ 
ments/s.  Its  signal-to- noise  ratio  exceeds  anything  available  commercially  or  other¬ 
wise  at  comparable  resolution  and  framing  rates.*  The  KFS  CCD  was  designed  in 
collaboration  with  JPL  (M.  Wadsworth,  S.  A.  Collins,  and  S.  T.  Elliott)  and  fab¬ 
ricated  by  Lockheed  Martin  (Milpitas,  Palo  Alto).  The  camera  head,  electronics, 
data-acquisition  system,  and  computer  interface  hardware,  firmware,  and  software 


8  Now  with  PixelVision,  Inc. 

*  The  conversion  dynamic  range  is  programmable,  with  13,  14,  15,  or  16  bits/pixel  at  lower 
framing  rates. 
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were  designed  and  fabricated  at  Caltech  by  Dan  Lang.  Other  KFS  system  design 
and  components  were  developed  in  collaboration  with  Pavel  Svitek  (Aeronautics), 
Steve  Kaye  (Physics),  Brian  Kern  (Astronomy),  and  Paul  Dimotakis. 

In  its  present  implementation  (Phase  1),  the  system  is  capable  of  160  fps,  corre¬ 
sponding  to  1.6  x  10s,  12-bit  measurements/s.  At  this  rate,  it  can  record  4000 
to  12,000  frames,  depending  on  data  signal-to-noise  ratio  by  utilizing  a  front-end, 
lossless  data  compression,  before  offloading  to  hard-disk  storage.  At  25  —  50  fps,  it 
can  record  for  minutes.  Upgrades  in  progress  will  increase  the  rate  to  its  1000  fps 
(Phase  2)  and,  subsequently,  its  recording  capacity  at  the  highest  rate  to  millions 
of  frames  (Phase  3). 

It  was  preliminarily  used  (prior  to  Phase  1  completion)  on  the  200"  Palomar 
telescope  (May  2000)  for  atmospheric-turbulence/aerooptics  measurements  using 
a  bright  star  as  a  reference  (with  Prof.  Chris  Martin  and  Brian  Kern;  Caltech 
Astronomy/Physics).  At  present,  further  increases  in  high-framing  rate  LIF  appli¬ 
cations  would  be  limited  by  available  high-energy,  high-quality-beam,  pulsed  lasers. 
Long-term  experimental  goals  include  3-D  unsteady  measurements  in  turbulence 
and  mixing,  combustion,  aerooptics,  and  many  other  areas.  This  is  a  vital  tech¬ 
nology  development,  broadly  useful  to  many  other  areas,  both  within  and  outside 
AFOSR  and  DoD  interests. 


1.7  Conclusions 

The  research  in  this  program  is  comprised  of  several  inter-related  parts. 

The  numerical-simulation  (Sec.  1.1.1)  and  experimental  studies  (Sec.  1.1.2)  of  hydro¬ 
carbon  flames  are  motivated  by  the  need  to  characterize  hydrocarbon  flames  un¬ 
der  variable  strain-rate  conditions.  Such  conditions  are  encountered  in  anticipated 
SCRAMJET-propulsor  configurations  and  complex-geometry  internal  flows,  which 
are  the  object  of  the  studies  outlined  in  Sec.  1.3.  In  SCRAM JETs,  in  turn,  many 
fuel-injection  scenaria  anticipate  transverse-jet  configurations,  which  were  studied 
as  outlined  in  Sec.  1.2.  Injected  fuel  can  be  in  the  form  of  droplets;  sphere/droplet 
wakes  in  the  relevant  Reynolds-number  range,  the  first  part  of  a  study  of  which  is 
in  ??,  which  continues  under  our  current  AFOSR  grant  F49620-01-1-0006.  Flow 
with  density  inhomogeneities  that  is  subject  to  acceleration  forces  in  the  direction  of 
density  gradients  —  virtually  assured  in  inhomogeneous  turbulent  flow  in  both  in¬ 
ternal  and  rotating-machinery  flows  —  is  subject  to  the  Rayleigh-Taylor  instability. 
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This  has  been  the  subject  of  investigations  in  other  communities,  but  not  by  much 
of  the  fluid-dynamics  community.  It  is  the  subject  of  direct-numerical  simulation 
and  a  preliminary  experimental  effort,  as  outlined  in  Sec.  1.4. 

Our  modeling  of  hydrocarbon  flames  and  combustion  led  to  the  computational  and 
experimental  work  on  low-C  hydrocarbon  flames  (Sec.  1.1)  and  highlighted  signifi¬ 
cant  deficiencies  in  the  ability  to  simulate  such  phenomena.  Low-C  hydrocarbons 
are  the  product  of  thermal  cracking,  either  in  flames  or  as  anticipated  in  the  en¬ 
dothermic  use  of  the  fuel  in  hypersonic  vehicles.  Their  relation  to  ignition,  flame 
stability,  and  extinction  are  important  to  air-breathing  propulsion,  in  particular, 
and  the  AFOSR  mission,  in  general.  Under  continuing  AFOSR  support,  we  are 
presently  engaged  in  a  program  to  address  these  deficiencies  as  well  as  provide  a 
substantial  increase  in  the  available  experimental  database  to  assist  in  the  improve¬ 
ment  and  validation  of  chemical-kinetics  models.  A  first  presentation  (Vagelopoulos 
&  Dimotakis  2001)  of  our  results  will  take  place  at  the  up-coming  2nd  U.S.  Joint 
Meeting  (25-28  March  2001,  Oakland,  CA). 

At  this  time,  much,  if  not  most,  of  turbulence  theory  relies  on  assumptions  of 
(at  least,  local)  isotropy  and  homogeneity  to  reduce  the  number  of  variables  that 
must  be  modeled.  Our  work  on  transverse  jets  (Sec.  1.2),  an  important  flow  in  its 
own  right  and  in  many  high-speed  air-breathing  propulsion  and  turbulent-mixing 
contexts,  as  noted  above,  highlighted  the  need  for  descriptions  of  turbulence  that 
do  not  rely  on  these  assumptions. 

The  same  conclusion  was  arrived  at  on  the  basis  of  studies  of  Rayleigh- Taylor- 
instability  flow  (Sec.  1.4),  as  well  as  a  host  of  others  on  mixing.  The  latter  illus¬ 
trated  significant  similarities  with  mixing  in  high- Re  shear  layers  as  well  as  with 
turbulent  microscale  behavior  of  turbulent  jets.  One  of  the  largest  direct-numerical 
simulations  of  this  phenomenon  was  conducted  as  part  of  this  study,  relying  on  the 
full  Navier-Stokes  equations,  including  molecular  diffusion  between  the  two  partic¬ 
ipating  fluids.  It  was  undertaken  in  collaboration  with  the  Lawrence  Livermore 
National  Laboratory  and  performed  on  DOE’s  ASCI  computers. 

Advances  were  also  realized  in  flow-field  velocimetry  (Sec.  1.5),  based  on  previous 
work  on  Image  Correlation  Velocimetry,  as  well  as  on  high-performance  digital¬ 
imaging  technology  (Sec.  1.6).  These  developments  are  crucial  to  the  success  of 
several  AFOSR-supported  research  projects  at  this  time,  under  AFOSR  grants 
F 49620-00-1-0036  and  F49620-01-1-0006. 
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The  research  program  undertaken  as  part  of  this  grant  was  very  productive.  A 
large  number  of  publications,  reports,  presentations,  interactions  and  collaborations 
stemmed  from  it,  as  outlined  in  Secs.  4  and  5. 
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Progress  and  events  not  documented  in  the  most-recent  annual  report  (31  August 
2000)  that  occurred  prior  to  the  end  of  the  present  Grant  (30  November  2000) 
include: 

-  29  August  2000:  P.  E.  Dimotakis,  “Recent  advances  in  turbulent  mixing,”  20th 
IUTAM  Congress ,  22  August  -  2  September  2000  (Chicago,  IL).  Invited  sec¬ 
tional  lecture  (Paper  EW1). 

-  11  September  2000:  P.  E.  Dimotakis,  “Progress  in  turbulent  mixing.”  Invited 
lecture,  Venderbilt  U. 

-  12  September  2000:  P.  E.  Dimotakis,  L.  Ribarov,  R.  Pitz,  S.  Nandula,  J. 
Wehrmeyer.  Discussions  on  further  collaboration  on  applications  of  Image 
Correlation  Velocimetry  to  gas-phase,  molecular-tagging  techniques. 

-  3  October  2000:  J.  Bergthorson,  “An  experimental  investigation  of  C1-C2  hy¬ 
drocarbon  flames,”  GALCIT  Fluid  Mechanics  Research  Conference  presenta¬ 
tion. 

-  3  October  2000:  M.  Wang,  “Statistical  and  Spectral  Analysis  of  High  Schmidt 
Number  Scalars  in  Grid  Turbulence,”  GALCIT  Fluid  Mechanics  Research  Con¬ 
ference  presentation. 

-  20  October  2000:  P.  E.  Dimotakis,  “Fluid  and  Solid  Mechanics,”  Workshop  on 
Barriers  to  Predictive  Simulation  in  Science  &  Engineering,  Lawrence  Liver¬ 
more  National  Laboratory. 

-  3  November  2000:  P.  E.  Dimotakis,  “Turbulent  drag  reduction,”  DARPA  out- 
brief  (Jason  study). 

-  7  November  2000:  M.  B.  Johnson,  “Mixing  and  Combustion  in  Complex  High 
Speed  Flows,”  GALCIT  Fluid  Mechanics  Research  Conference  presentation. 

-  19  November  2000:  J.  W.  Shan  and  P.  E.  Dimotakis,  “Turbulent  scalar  mixing 
in  transverse  jets,”  Am.  Phys.  Soc.,  53nd  Annual  Meeting,  Division  of  Fluid 
Dynamics ,  (19-21  November  2000,  Washington,  DC),  BB.007. 
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Consultative  and  advisory  functions  to  other  labs,  government  agencies,  especially 
AF  and  DoD  (include  institutions,  locations,  dates,  names): 

-  Dimotakis,  P.  E.:  Lawrence  Livermore  National  Laboratories.  Consulting  on 
compressible  turbulence  and  inertial-confinement  fusion  (1993  to  present). 

-  Dimotakis,  P.  E.:  JASON  group,  MITRE  (1997  to  present).  Consulting  to 
DARPA  and  other  government  agencies. 

Transitions  (cases  where  knowledge  resulting  from  this  research  has,  is,  or  will  be 
used  in  a  technology  application): 

-  Contributions  to  High-Energy  Laser  (HEL)  development  (1970-90).  Contact: 
Dr.  Wilhelm  Behrens  (310  812-0468),  TRW  Rl-1044,  1  Space  Park,  Redondo 
Beach,  CA  90278. 

i.  Shear-layer  growth  and  entrainment  data/models. 

ii.  Reynolds-number  effects  on  mixing. 

iii.  Use  of  gas-mixtures  for  index-or-refraction  matching. 

iv.  Laser-induced  fluorescence  diagnostics. 

New  Discoveries,  inventions,  patents:  None 


6.  Honors  /  awards 

Honors,  degrees,  and  awards  received  during  period,  as  well  as,  ‘lifetime  achievement 
honors  such  as  Nobel  prize,  honorary  doctorates,  and  society  fellowships  prior  to 
this  effort: 

•  Bond,  C.  L.: 

-  Ph.D.,  California  Institute  of  Technology  (June  1998). 

•  Dimotakis,  P.  E.: 

-  John  K.  Northrop  Chair,  Aeronautics,  Caltech  (February  1995). 

-  Associate  Fellow,  AIAA  (June  1989). 

-  Fellow,  Am.  Phys.  Society  (November  1980). 
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•  Leonard,  A.: 

-  Fellow,  Am.  Phys.  Society. 

-  Vice-Chair,  Am.  Phys.  Society,  Division  of  Fluid  Dynamics  (1997). 

-  Chair,  Am.  Phys.  Society,  Division  of  Fluid  Dynamics  (1998). 

-  Theodore  von  Karman  Chair,  Aeronautics,  Caltech  (December  1999). 

•  Slessor,  M.  D.: 

-  Ph.D.,  California  Institute  of  Technology  (June  1998). 


